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Abstract 

We describe the Kerr black hole in the ingoing and outgoing Kerr-Schild horizon penetrating 
coordinates. Starting from the null vector naturally defined in these coordinates, we construct the 
null tetrad for each case, as well as the corresponding geometrical quantities allowing us to explicitly 
derive the field equations for the VPo^ and perturbed scalar projections of the Weyl tensor, in- 

cluding arbitrary source terms. This perturbative description, including arbitrary sources, described 
in horizon penetrating coordinates is desirable in several lines of research on black holes, and con- 
tributes to the implementation of a formalism aimed to study the evolution of the space time in the 
region where two black holes are close. 
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1 Introduction 



A particle falling unto a black hole is an efficient mechanism for generation of gravitational waves. This 
problem is not only of great importance to gravitatonal wave astrophysics, but it also has important 
application of numerical general relativity. 

The fact that several gravitational observatories (LIGO, VIRGO, GEO600, and later LISA) will start 
to be working soon is one of the most exciting news lately received by the general relativistic community. 
Those observatories will allow a much deeper insight on the phenomena near the black holes such as the 
above mentioned particle falling unto it, and later, will allow us to see beyond the recombination era, 
plus all the not yet expected ones. 

These technological achievements, together with a strong development in the computational tech- 
niques, definitely boost the great theoretical knowledge acquired in the last century. In particular, the 
theoretical study on gravitational perturbations goes back to Einstein and has produce excellent works 
such as the one by Chadrasekhar M or the one by Frolov for just to mention a couple of them. Also, 
there is the work by Teukolsky || who, using the Newman Penrose formalism for null tetrads (see @, 
was able to obtain an equation for the ingoing and outgoing radiative parts of the perturbed Weyl tensor, 
namely 'J'o' 1 ' and 'J 4 , respectively, which were separable in the frequency domain. The radial part of 
those equation were useful in obtaining information about the asymptotic behavior of the gravitational 
perturbations in the vacuum case. 

The complexity of the equations prevented further advances until the development of numerical tech- 
niques for studding the evolution of the equations. There was large effort in understanding the collision 
of black holes evolving directly the Einstein equations 0, 9, as well as those studies comprising semi 
analytical results where they used the perturbation equations derived by Teukolsky with two important 
changes: First, with a coordinate system which avoids the coordinate singularity, and second, they worked 
in the time domain, to be able to study the evolution. 

Campanelli et al. |(| presented a work in this last direction and obtained the evolution of the VE^ 1 -* 
Weyl tensor for the case without sources. In the present work, we continue within such line of research, 
obtaining the equations which allow the study by numerical evolution of the and \I/4^ perturbed 

Weyl scalars in the case with arbitrary sources generating the perturbation. Furthermore, families of 
initial data with appealing features such as describing regularly multiple moving spinning black holes J?J , 
use the Kerr Schild horizon penetrating coordinates. 

We stress the fact that in the present work we explicitly derive the field equations for the and 
^4^ perturbed scalar projections of the Weyl tensor, including arbitrary source terms. Even though 
there is the coordinate transformation from the Boyer Lindquist to the Kerr Schild ones, the change of 
tetrad involves rotations and dilatation of the null tetrad which prevent to treat the Teukolsky equation 
as a tensor operator, and it is necessary to derive the equations from their general definitions, as is done 
in the present work. 

The work is composed as follows. In the next section we derive the Kerr metric in ingoing and outgoing 
Kerr Schild coordinates, and present the field equations for and \EV , deriving the master field 

equation. Next, we take the naturally defined null vector and construct the null tetrad, the differential 
operators, the spinor coefficients and the Weyl scalars and, describing the main steps in the deduction, 
we finally derive explicitly the expression for the perturbation equations for arbitrary sources valid in a 
Kerr Schild description, in the ingoing as well as in the outgoing formulation. We test our result in simple 
examples and conclude mentioning several case where will be of interest apply the obtained equations. 

2 Kerr Schild line element and master perturbation equation 

The rotating black hole is usually described by the line element in Boyer Lindquist coordinates || : 

ds 2 = A ^ _ as . n2 _ sin^ ^ 2 + ^ ^ _ adf y_ E /d£ + dQ2 \ ^ (1) 

V AA / 
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2Mr + a 2 , and £ = r 2 + a 2 cos 2 6. 
However, as mentioned above, this description has a divergence at the external horizon: r + = M(l 



+ (-p) 2 ) , which prevent that the evolution analysis could be performed in the vicinity of the external 
horizon. Moreover, it has been established that the gravitational wave emission of several phenomena 
such as the matter accretion or the inspariling of an object unto the black hole, occurs precisely within 
this region (9). Detailed studies show that there is a sweep upward in frequency, during the last 15 
minutes before the collision and final coalescence, that is, precisely in the region close to the external 
horizon for which the Boyer Lindquist coordinates are not well suited. 

Thus it is imperative to be able to perform the description of the phenomena within a coordinate 
system which does not have those coordinate problems. Whit in this spirit the Kerr Shild description not 
only avoids the problem of coordinate singularities at the horizon but is the most natural description to 
be matched with the null tetrad formulation. In this way, following 0, we define the coordinates t and 
if in terms of the outgoing and ingoing ones, u in : 



r 2 + a 2 a 
du in = dt ± dr, dip in = dip ± -rdr, (2) 

out /\ out /\ 

and, defining in each case a new temporal coordinate 

dt in = du in =F dr, (3) 

out out 

We obtain that the Kerr space time is described as: 

j2 , 2Mr, , 2 4Mr, , 4Afrasin 2 6> , / 2Mr\ , 2 
ds 2 = (1 )dt 2 i„ =F — — — dt i» dr H dt u> dtp m - H dr 2 

out Jj out J] out out \ y] J 

( 2Mr\ , A , , ,«■>[■> 9 2Mra 2 sin 2 6»\ , „ , , 
±2 1 + J a sin 2 Odrdtp u, - Sd6» 2 - ( r 2 + a 2 + J sin 2 6dtp 2 ,„ , (4) 

which is the known Kerr Schild description. Noticing that this last expression for the Kerr line element, 
Eq. (||) can be rewritten as: 

2Mr 

ds 2 = (rj/a, in — Ifj, i„ l v in )dx^ dx" , (5) 

out 2_j out out 

with ri U u s« is the flat space metric, which in this case takes the form: 

out 

r/nv in dx^dx u = dt 2 i„ - dr 2 - Hd9 2 - (r 2 + a 2 ) sin 2 6dip 2 i„ ± 2 a sin 2 Odrdtp ,. , (6) 

out out out out 

and / M in is a null vector, null as much for the flat metric as for the complete one, and is given by: 

out 

= (1, ±1,0, -asm 2 0). (7) 

We want to stress the fact that, as mentioned above, the line element in the Kerr Schild coordinates, 
Eq.(^), not only has no geometric divergence all the way from infinity to the interior of the black hole, 
but also defines naturally one of the null vectors. 

Thus, we use this null vector, defined by Eq.(Q), to construct the null tetrad, the differential operators, 
the spinor coefficients and the Weyl scalars, which are needed to study the gravitational perturbations 
within the Newman Penrose formulation. The null tetrad l,n,m and m* is defined by the following 
identities: 

I ■ n = 1 and m ■ m* = — 1, (8) 
and from the relation of the tetrad with the metric coefficients: 

= 2[Z( M n„) - m^m* u )]. (9) 
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Before presenting explicitly the ingoing and outgoing cases, we derive the field equations for the 
perturbed Weyl components, valid for any algebraically special type D space time. 

Perturbing the Bianchi identities in the NP formulation ||, Hi; taking into account the expressions 
derived in the above section, coupled equations for the perturbed Weyl scalars 'J'c/ 1 ' and Vl/i^ and 
$4^\ are obtained. Using the generalized commutation relations between the spinor operators given in 
the appendix A, Eq.(|37|), as well as the perturbed Ricci identities obtained form jj], ^| and the expressions 
derived in the above section, the following separate equations for the perturbed Weyl scalars ^o' 1 ' an d 
^4^-* are obtain: 

{(D - 4p - p* - 3e + e*)(A + p - 4 7 ) -{8-2(3- At){5* + 4/3* - 3vr) - 3* 2 }*o (1) 

= 87r 2 T^T^W, (10) 



{(A + 4p + p* + 3 7 - 7 *)(D - p + As) - (S* - 2/3* + 7tt - t*)(6 + 4/3 - t) - 3f 2 } 



= 87r_ 2 TT |u/ W, (11) 
where S T F are given by: 

aT"" = -{5 - 2/3 - 4t)(<5 - 7T*) l v - (D - 4p - p* - 3e + e*)(D - p* - 2e + 2e*) m" m u 
+ [(S - 2/3 - 4t)(D - 2p* - 2e) + (D — 4p — p* — 3e + e*)(S - 2/3 + 2tt*)] & m v \ 
-aT"" = 

-{5* - 2/3* + 7tt - t*)(5* +2tt- t*) n* 1 n v - (A + 4/i + p* + 37 - 7 *)(A + p* + 2 7 - 2 7 *) to* m to*" 
+ [(<5* - 2/3* + 7tt- t*)(A + 2p* + 2 7 ) + (A + Ap + p* + 3 7 - 7 *)(<5* - 2/3* + 2tt - 2r*)] m* v) , 

(12) 

with T^ 1 ' is the source of the perturbation. These equations were first derived by Teukolsky in 1972, 
H and, as long as he was interested in the asymptotic behavior, the Eddington Finkelstein coordinates 
which he used, were well suited for the analysis. 

The studies are performed on the \E f o ( ' 1 ' ) an d 4 , 4^ 1 ^ perturbed Weyl scalars as long as those are the 
most significant components in the gravitational radiation || . 

Acting on p -4 ^f^ 1 ' and using the commutation relations between the operators, Eq. (p7[) , and the 
action of the operators on the spinor coefficients, Eq.(p9|) into the appendix A, the perturbation equations 
for * (1) and p- 4 * 4 (1) , Eqs.(|ll|), can be written as a single master equation in terms of the parameter 
s: 

{AD - 5*5 + {p* - {2s + 1) 7 - 7*) D - ((2s + l)p + 2se) A - (2(s + l)/3* - (2s + 1> - r*) 6 
+ (2s/3 + (2s + 1)t) S* + s (2s + l)[2(p + e) 7 - 2/3 tt - 2ttt + 2/3* r - * 2 ] - 

2 s [(Ae) - {5*/3) +e{p*- 7*) + /3r*] + 4s (s + 1) /3/3*} s i\> = 8tt p^ S T» V (13) 

with 2ip = y -2ip = P~ 4 ^4^- Even though the ^o^ 1 ' and ^4^ perturbed Weyl scalars are 

the most significant components in the gravitational radiation, it would be interesting to analyze the 
evolution of the rest of the perturbed Weyl scalars, ^i*- 1 ^, ^2 and "fa*- 1 -*. Some studies indicate that 
it is possible to obtain a single master equation for the five of the perturbed Weyl scalars |jlC 
These studies are currently under investigation and will be reported elsewhere. 
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3 Null tetrad formulation and explicit derivation of the pertur- 
bation equations for V 1} and ^ 4 (1) with sources 

Next, we presents the main steps in the derivation of the explicit form of the field equations for the 
perturbed and ^4^' Weyl components with sources. In the appendix B, we present the main 

expressions for the action of the operators on the spinor coefficients, needed in the derivation of the 
explicit equation. We are going to present both cases, the ingoing and the outgoing. 
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3.1 Ingoing case 

Starting form the naturally defined null vector, given by Eq.(|7|), using the equations Eqs.^J^), we con- 
struct the following ingoing null tetrad: 

,0, -a 

(A, r 2 + a 2 (cos 2 6 - sin 2 0) + 2Mr, 0, -aA sm 



i Min = (f, 1,0, -a sin 2 i 
Mm 2S 

m u; „ = — — — [— i a sin 0, — i a sin#, E, i (r 2 + a 2 ) sin 0] . (14) 



V2 (r + io cos 
From which we obtain the covariant null vectors: 

fto = (1,-1,0,0), 



n- 



= — (r 2 + a 2 + 2Mr,A,0,2a) 



V2 (r + ia cos 0) 
as well as the corresponding differential operators: 



2E 

= (iasin0,O,l,icsc0), (15) 



D in 


= ^hA = d tiD - d r , 






A in 


= = ^ [(r 2 +a 2 + 2Mr)d ti 


1 + Ad r 4 


- 2ad ip 


Sin 


= m^indu = — ^(ia sin 9d t , + de + 
V2 


% esc 9d Vi 


J, 



with the spinor coefficient p in defined below in Eq. (|l~7|) . 

The spinor coefficients are computed directly from their definition ||, |^] , obtaining: 

1 , A (r-M) 

2ej n pi n — p in , /i; n — — pi n , 2^1!, — p- ln + /i ; n — . 



(16) 



r + ta cos 2 E 

_ basing _ cotg „ _ _ _ msinfl 2 

Tin y^2S ' ' \Z2 ^ n n — ^~* n n ' ^ n — \Z2 ' 

and finally, the only non zero Weyl scalar is: 

* 2in = -Mp-J. (18) 
Using the explicit form of the differential operators, Eqs.(fLq), of the spinor coefficients and the Weyl 



scalar Eqs. fll7| , |18|), for the null tetrad Eq.(|15|) of the Kerr metric described in Kerr Schild coordinates, 
Eq. (|J) , with the help of the expressions for the action of the operators on the spinor coefficients given 
in the appendix A, Eq. (|39|) , (|4C|), we derive the equations for the perturbations of the ^o' 1 ' and "J^ 1 -* 
perturbed Weyl scalars due to arbitrary source terms in ingoing coordinates: 

d 2 d 2 d 2 d 2 d 2 1 d 2 



dt 2 n drdt in dr 2 ~ drdip in d9 2 sin 2 6 dipf n 

2 ((« + l)(r - M) + 2sAe in ) — - 2(sr + (s+ 1)M - se- m (E +4Mr)) — 

(Jf ^'t'in 

/ « „ « * d n /.cos6> \ S 2rt . 16a 2 .„ 2Mrcos 2 6l 

- cot0 + 4srtan6le ill — -2s i — _ + 4ae in s(l - scot 2 0) + — -(1 

06 \ sm J oipi n E E 

-4*6* ^( S + l)r+(r 2 ~s{r 2 -a 2 cos 2 9))^]^^ l&nll^T^, (19) 
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where s T' lly are given by the following expressions: 



2 T^ = 2 T u l*l v + 2 T mm rn>*m v + 2 T lm Fm v 
_ 2 T^ = _ 2 T nn n' i n 1 ' + - 2 T m * m *m*> J 'm t- _ 2 



+ _ 2 T nm n^m*" 



(20) 



with 



n 2 r) 2 
2i cos 9 9 



2ia sin# 



d9dt ir 



2a 



if 



dtp in dt in 



9 2 2 
dip? n dd 2 
9 9 . 



2i esc 9 



d9dip- u 



sin 2 6 dipi n 
ia sin 8 p* n }}, 



.9 „ . .6r — 2iacos9 , , . 9 9 n 9 

c °t 9— + 2ia sm9 — ha sin 9— h — + i esc 9— — 

09 E OT in 09 difu 



dr 



6r — I0iacos9 9 
E )( 9i~ 



9 4 

— ) H ^(r 2 — 5a 2 cos 2 — 4riacos#)}, 

ar E^ 



-2pin rr . „ d 9 

7T {(wsm t + i 



,9 .„ iasin#. , 9 
>- h cot# + (2r - 3zacos6») )(-— 

dipin E Stir. 



— ) 

dr' 



dp- ln (ia sinf? 



9 9 9 9 

c °t#— + 2iacos9— + iap* n sin9[iasm " 



9 d_ 
dt~ + dl 

9* 
dtl 
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difu 
d 2 
89dU n 



cot ( 



2a- 



ia sin 0(0;, 



(6r 2 + 5iar cos 6> - 5a 2 cos 2 9) } , 



d 2 



t) 2 



CSC 



9 2 „ „ d 2 

W + 2lcsce d9bV v 
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89 



-^{A + (A + 



at 



dt m 



2A 



d_ 

dr 



4a- 
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dt m 
d 



d9 h 



l CSC ( 



difin dt h 



dipin 

a / a 5 , 9 . d 
A(A— +4a-— ) — 

Or 0<Pin <JT 



»> Pin » 



2a 2 sin 2 6»p? nJ , ul 



4a^ 



4M(3r 2 + a 2 )-^- - 2(r - 3ia cos 0) | (A+ J- + a|- + 2a-^-) 
OT in E 9t in 9r 9<^ in 



<itlCOS ypin-jrlPin 4 , 



A 2 
E 



Pfn r / a 9 . 9 „ 9 



A 9 9 

(2r + dia cos 9) r + M)(iasm.6— — 

E 9<in 9(9 



. 9 . , „ ,„ . „. lasmO... 9 .9 „ 9 . „ . „ A 
icsc9— ) + (cots' — (dr + iacos9) )(A + — h A— + 2a— J + oiarsm9p in — 

Oipi-n E 9t; n Or 0<Pm ^ 



2(r - M)(cot0- (dr + ia cos 9) mS ™ 6 ) - (r + diacos9) cot 9^}p hl 4 , 



(21) 



where we have also act the ^T* 11 ' operator on p . 

Once obtained these perturbation equations for arbitrary sources within a description free of geometric 
divergencies, it is possible to proceed further and perform the numerical evolution of the and 
perturbed Weyl scalars for cases as interesting as the ones due to an orbiting object, or the one mention 
above of an inspiring and finally colliding and coalescing object, or those due to accreting matter, for 
mention some examples. These studies are currently under way. 

We finish this presentation testing the perturbation equations Eqs.(|l9|) in Schwarzschild described 
with the Kerr Schild coordinates, which implies setting the angular parameter a — 0, and the Eqs.( p^| , pT| ) 
reduce to the form: 



9 2 9 2 



drdU 
9 
dls 



-A 



9 2 



d 2 



1 



dr 2 89 2 sin 2 9 dip. 
9 . cos 9 9 



^i-2(( S +l)(r-M) + 2 S A Ein )|- 



-2 (a r + (s + l)M) - cot 9 ^ - 2si 7" 2 " " - s(l - scot 2 8)] s ip = 16ttE s T» v T m „ (1) , 
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(22) 



with A + = r 2 + 2 M r, the angular coordinate ip remains unchanged in this case, and now A = r 2 — 2 M 1 
The operators S P" are given by Eq.(pb]) but now the projected operators reduce to: 



1 r 2 „ d 2 d 2 „ n d 2 2icos6 d J, 
2r 2 dipf n d9 2 96d(p m sir 6 oip in dO J 

T mm -U— - — --)(— -—)-—\ 

X{ dt in dr r )[ dt in dr' r 2 h 

r™ 2 r . 8 d 3 . , d „ d 

* = -^ {( ^--^~^ )( ^ +lCSC ^ +COt ' )} ' 

"* m * = -^{A+(A+— + 2A— )— + A 2 — + 12Mr 2 — -2-(A+— + A — )}r 4 . 

4r 4t +l + 9t in dr'dt ia dr 2 dt m r l + di in dr" ' 

— ■• 1 r/^A 9 9 a \ / d „ a 

r{ +A-+2--r + M -- 



\/2> 3 L V di in 0r r / V <96> dip- 

( 8 8 \ A 
+ cot6» A+— + A— - 2(r — M) cot 9 — cot 6 — }r 4 . (23) 
\ at in ar / r 

We try the simple case when the perturbation is due to dust radially falling unto the black hole, for 
which the stress energy tensor is given by 

T" v = d(r, ', (24) 

with d(r, t) the density and the four velocity, which for the radially falling case has only radial and 
temporal components: 

u" = (u*,u r ,0,0). (25) 

Projecting the stress energy tensor, Eq.(f24|), unto the null tetrad which for the Schawrzschild case is 
given by: (from Eq.(|l5|) with a = 0) 

i„ = (1,-1,0,0), 

n„ = ^(A,-(r 2 + 2Mr),0,0), 

m p = (0,0, l,t sin 6»), (26) 

V2 



and taking into account the components different from zero of the four velocity, Eq.(|25|), we obtain that 
the stress energy tensor projections unto the null tetrad are, 

T„ = «*»(»•, i)(«°-u r ) a . 

T nn = ^M(Au°-A + ^) 2 , 

Tim Tmm T n . m * -^m*m* 0. (27) 

Notice that the non zero projections are functions of r and t only. Thus, substituting the projections in 
the rhs of Eqs.(p2|), there is only left the action of the operators on Tu and on T nn . But those operators 
for the Schwarzschild case are purely angular operators, and as the projections are functions of (t, r), the 
action is zero. In this way, we conclude that radially infalling dust into the Schwarzschild black hole does 
not produce perturbations on the and perturbed Weyl scalars for both cases. 
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This result is expected considering that the gravitational waves carry quadrupolar radiation and 
higher, and that the radially infalling dust into the Schwarzschild black hole will produce only monopolar 
perturbations, which will not be reflected in the and \&4^ perturbed Weyl scalars. 

It can be seen that the same happens with a radially infalling massive scalar field, with stress energy 
tensor given by 

T„ v = 4>» <t> v - \g»u {4> a <t> a - 2V (<!>)) , (28) 

with the scalar field a function of (t, r) only, and where subindex and superindex imply derivative with 
respect to x M . Again, the only non null projections needed in the perturbation equation are 

T„ = (0 + 0') 2 and T nn = ^(A+ - A 0') 2 , (29) 

where dot and prime stand for derivatives with respect to t and r respectively. Once more, they are 
only functions of (i, r), thus the action of the operators on them is zero, so the massive scalar field 
radially infalling unto the Schwarzschild black hole does not produce perturbations on the 'Jo*' 1 '' and 
*S>4^ perturbed Weyl scalars. 



3.2 Outgoing case 

In this section we present the analogous derivation of the field equations for the perturbed scalar pro- 
jections and ^E , 4 ( - 1 ' ) of the Weyl tensor, including sources. The procedure, although length, is 
straightforward and makes the corresponding steps as in the ingoing case. Thus, from the null vector 
naturally defined in the Kerr Schild form, Eq. ([?]), we construct the following contravariant null tetrad: 

^out - (1, -1,0, -asm 2 0), 

n Mout = -^(A,-(r 2 + a 2 (cos 2 0-sin 2 0) + 2Afr),O,-aAsin 2 0), 

77 Vout = ~~r [— i a sin i a sin#, £, i (r 2 + a 2 ) sin#] . (30) 

i U V2 (r + i a cos 8) 

Thus, the covariant null tetrad is: 

P out = (1,1,0,0), 

n^out = ^(r 2 + a 2 + 2Afr,-A,0,2a), 

(ias'mO, 0, 1, icsc9). (31) 



V 2 (r + i a cos 

From which we obtain the corresponding differential operators: 
D ou t = I out^i = c*t out + d r , 

A out = ^^8^ = ^[(r 2 + a 2 + 2Mr)d tout - A8 r + 2ad Vout ] , 

Sou* = m^fi^ = ^ (to sin 9d tout +8g + i csc 6d VoM ), (32) 
V2 



with the spinor coefficient p ou t defined below in Eq.(|33|). 

The spinor coefficients are computed directly from their definition |^] , obtaining: 

K out — Cout = -^out = ^out = £out = 0, 

1 A r-M 



Pout . „ , Mout n „ Pout , Tout A* ' n ! 

r — i a cos 2L 2 2j 

_ ia sinO 2 _ cotfl „ _ _ _ msinfl 

T^out ~~ ~^2, ^ out ' Pout — 2\/2 ^° U * ' ^ out — ^out P out J Tout ~~ \f% S ' 
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and finally, the only non zero Weyl scalar is: 

*2out = Mpout 3 - (34) 

As mentioned above, the Teukolsky equation, that is the field equation the perturbed scalar projections 
'Jo'' 1 '' and of the Weyl tensor, including sources, Eq.(|l3|), as well as the expression for the operators 

action on the sources, Eq.(|2C|), are valid also for the outgoing case. Thus, using the geometric quantities 

39 4l]) after a quite lengthily 

algebraic manipulation, we obtain the following equation for and $4 ^'perturbed Weyl scalars due 

to arbitrary source terms in outgoing coordinates: 



given by Eqs.(B2 133,34), and their properties given in the appendix B, Eq.( 



[(S + 2Mr)X 



4 Mr 



drdt 



8 2 

dr 2 drdip ou t 



o 2 



d 2 



dO 2 sin 2 9 dip, 



out 



-2(s + l)(r — M)— + 2(sr + (s + l)M + s ia cos ( 



9 cQ 9 
dt out 09 



x cos 6 8 
-2(s + l)i — 5-^ — 
sin 9 dip ox 



s(i - s cot 2 e)] a i> = i6tt e s t^ 



(1) 



(35) 



where are given again by Eq.([20|), but using the outgoing terms, the projected operators have now 

the following expressions: 



rrll 

2± 



Tim 



*2 
P out 



9 9 8 8 
ia sin 6 (ia sin 6— h 2— + 2i esc 9— )— — 

OTout 09 Olfiout Otou 



Sin ^ (CSC ^ } 



8 8 d 2 

+2i — (csc9- )-csc 2 9— 

d9 otfout dip. 



2 

out 



iasinf? . 8 

2 (Sr + iacos9 wsmff- — 



9_ 

"09 

-( 



9 



S 



9iout 
P out 



[(■ 



9<^out 

9 3r + ia cos 
1 

dr 
d 



ia sin p* ut )], 
9 



-)(" 



V2 9i ou t 
— 2zasm0(2p out 



S Stout 

9 3r + iacos0 
97 + £ 

5 . . 2 „ 



— ) 

dr' 

)(iasin( 



4 

£' 
9 



9i n 



_9 

96» 



l CSC ( 



9(p 



cot 9) 



8i, 



out 



p 2 

[ia sin 0(ia sin t/- 

2 OT ou t 



P ouJJ> 

8 



2^ 
96> 



2i esc 9 



d 







difout dt 



I) 2 



CSC 



iasin0 ._ _ . .. . n 8 

— — — (7r 4- Ilia cos 0)(zasm0 



d d 

-2i — (csc9 ) + 2p ut 

89 dfout 



a 2 sin 2 1 



9i ut 
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Again we recall that these equations were directly obtained working from the original field perturbation 



4 Conclusions 

Scratching form the top, we have derived the equations for two of the radiative parts of the perturbed 
Weyl tensor, namely v E'o ( ' 1 ^ an d in the Kerr spacetime described in the Kerr Schild coordinates, 

including sources. We have done son for the ingoing as well as for the outgoing descriptions. 

Even thought the procedure followed in this work is straight forward, we consider that the equations 
obtained are a needed step in the research line of semi analytical studies of perturbed black holes, and 
will be useful for community. 

The present description allows us to perform studies as near to the horizon as is wished, even inside 
the horizon for that matter, as long as the coordinates are regular. Also, there are several works which 
describe initial data, namely in the Kerr Schild description, so our work fits well with them. We remark 
the fact that the derived equations with sources presented in this work were obtained directly from the 
Bianchi identities within the general null tetrad formulation. The reason for this is that the spinors 
operators are not scalar quantities and the changes from one tetrad to another are not due to coordinate 
transformations only. Thus, it is not safe to use the equations for the perturbed components derived by 
S. Teukolsky, back in 1972 ||, where he used the Eddington Finkelstein coordinates (a description well 
suited for the asymptotic analysis, which is the one he performed), and use the coordinate transformation 
to describe those equation in the Kerr Schild coordinates. 

Thus, we derived the equations and tested them in simple cases where we obtain the expected results, 
such as the fact that radial perturbations do not generate radiation. The development of realistic situ- 
ations is currently under way. This work continues the program started by Campanelli et al. || where 
they studied the gravitational perturbation within a description free of coordinate singularities, without 
sources. We have introduce another null tetrad and derived the equations with arbitrary sources. 

The analysis and the equations presented in this work are a complement for those performed purely 
numerically. Joint efforts are currently under way. As mentioned above, there are several phenomena 
that can now be studied from this perspective, mainly the most important is to continue in the study of 
gravitational waves from a compact star in circular, inspiral orbits around a massive spinning black hole, 
with applications to observation by LISA and LIGO that is, to obtain predictions which could be 
detected and tested by the gravitational observatories which will start to work in the next year. 
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6 Appendix A 

The following commutation relations between the operators are repeatedly used in the work, and are 
valid for any algebraically special type D spacetime (we used that a — tt — /?*): 



equation, Eqs.(|l9|). 



D A = AD - (7 + 7*) D - (e + e*) A + (r + tt*) 5 + (r* + tt) 5*, 
55* = 5*5 + Oi-/i*)D+ (p-p*) A + (2/3-n*)5- (2(3* -n)6*, 
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(D + qp- p* - (p+l)e + e*) (5 - p/3 + qr) - (6 - p/3 + qr) (D + qp - pe) = 0, 
(A - q (i + p* + (p + 1) 7 - 7 *) (S* ~ p/3* + (p - g)7r) 

- (S* - P f3* + (p-q+ 1)tt — r*) (A — qp + py) = 0, (37) 

where p and g are arbitrary real numbers. This expression are valid for the outgoing and ingoing case. 
Finally, the following identities among the spinor coefficients are also useful: 

(3*T = —(3lt, 1T*p — Tp*, 

p*p = pp*, 2pj - 4/3ir - 2itt - * 2 = 0. 

(38) 

7 Appendix B 

In the derivation of the field equations, is useful to determine the action of the operators on the different 
spinor coefficients. Its explicit form is obtained from teh expressions for the spinor coefficients given by 
Eq. (|l7|, |33|) , and the operators given by Eq.(|l(], 32). Some of those actions are valid for both cases, the 
ingoing and the outgoing: 



Dp*=p* 2 , At* =-(/x + M*)t*, 5p = rp, 

S*p = -np, Sn* = 2((3 - n*) n* , S*p, = -(2n - t*) p , 

5p*=irp, 5*tt = 2(/3* - tt) tt, S*t* = (2/3* -it + t*)t* 



The following expressions are valid the ingoing case: 

Ve = \{p 2 - p- 2 ), D/3 = f(2/3 + r-r*), Dtt* = p(r* - tt), 

At* = -2/iV, Att = 7r*(/i + /i*), A/i = -0 + 2i)n, 

Ap* = -{p* + 2 1 )p\ A/3 = -^(2/3* +t-t*), S/x = M (2t - tt*), 

fe = f(r-T*), = (2t* — tt) p* , S*^- 1 = 2pr*€*. 



And, for the outgoing case, the following expressions are valid 

D/3 = /3p*, Dtt* = 2w*p*, At* = -{p + p*) t*, 

A7r = -2 7r^, A/i = -(//* + 2 7 ) M , A/i* = -(/i + 2 7 )/i*, 

A/3 = - M */3*, 5/x = p*{2ir* -t), <5>* = -(2tt* + r) p , 
5/3* = -p(3*. 



(39) 



(40) 



(41) 
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